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Abstract 

A new family of 2-component vector-valued coherent states for the quantum parti- 
cle motion in an infinite square well potential is presented. They allow a consistent 
quantization of the classical phase space and observables for a particle in this poten- 
tial. We then study the resulting position and (well-defined) momentum operators. 
We also consider their mean values in coherent states and their quantum dispersions. 
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1 Introduction 



Even though the quantum dynamics in an infinite square well potential rep- 
resents a rather unphysical limit situation, it is a familiar textbook problem 
and a simple tractable model for the confinement of a quantum particle. On 
the other hand this model has a serious drawback when it is analyzed in 
more detail. Namely, when one proceeds to a canonical standard quantiza- 
tion, the definition of a momentum operator with the usual form — ihd/dx 
has a doubtful meaning. This subject has been discussed in many places (see 
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[T] for instance), and the attempts of circumventing this anomaly range from 
self-adjoint extensions [lj to VT symmetry approaches [2]. 



First of all, the canonical quantization assumes the existence of a momen- 
tum operator (essentially) self-adjoint in L 2 (Kf) that respects some boundary 
conditions on the boundaries of the well. As has been shown, these conditions 
cannot be fullfilled by the usual derivative form of the momentum without the 
consequence of losing self-adjointness. Moreover there exists an uncountable 
set of self-adjoint extensions of such a derivative operator which makes truly 
delicate the question of a precise choice based on physical requirements [TH3] . 



When the classical particle is trapped in an infinite well of real interval A, the 
Hilbert space of quantum states is L 2 (A,dx) and the quantization problem 
becomes similar, to a certain extent, to the quantization of the motion on the 
circle S 1 . Notwithstanding the fact that boundary conditions are not periodic 
but impose instead that the wave functions in position representation vanish 
at the boundary, the momentum operator p for the motion in the infinite 
well should be the counterpart of the angular momentum operator L for the 
motion on the circle. Since the energy spectrum for the infinite square well 
is {n 2 , n G N*}, we should expect that the spectrum of p should be Z*, like 
the one for L without the null eigenvalue. This similarity between the two 
problems will be exploited in the present paper. We will adapt the coherent 
states (CS's) on the circle Pf5|6] to the present situation by constructing two- 
component vector CS's, in the spirit of [8], as infinite superpositions of spinors 
eigenvectors of p . 



In the present note, we first describe the CS quantization procedure. We re- 
call the construction of the CS's for the motion on the circle and the resulting 
quantization. We then revisit the infinite square well problem and propose 
a family of vector CS's suitable for the quantization of the related classical 
phase space. Note that various constructions of CS's for the infinite square 
well have been carried out, like the one in [9| or yet the one resting upon the 
dynamical SU(1, 1) symmetry [10j . Finally, we present the consequences of our 
choice after examining basic quantum observables, derived from this quanti- 
zation scheme, like position, energy, and a quantum version of the problem- 
atic momentum. In particular we focus on their mean values in CS's ("lower 
symbols") and quantum dispersions. As will be shown, the classical limit is 
recovered after choosing appropriate limit values of some parameters present 
in the expression of our CS's. 
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2 The approach via coherent state quantization 



Coherent state quantization |11|I12|I13|I14|I15II16|I17] is an alternative way of rep- 
resenting classical observables into a quantum system. The states used in it 
include Glauber and Perelomov CS's but lie in a wider definition that admits 
a large range of state families resolving the identity. Identity resolution is here 
the crucial condition. 

In fact, these states form a frame of reference well suited to represent classical 
quantities and, in that sense, work as a natural quantization procedure which 
is in one-to-one correspondence with the choice of the frame. The validity of a 
precise frame choice is asserted by comparing spectral characteristics of quan- 
tum observables / with data from the observational space. Unlike canonical 
quantization where the whole model rests upon a pair of conjugated variables 
within the Hamilton formalism [18], here we need the following elements. 

First of all let X = {x G X} be a set equipped with a measure /i(dx), and let 
L 2 (X, n) be the Hilbert space of square integrable functions f(x) on X: 

\f\ 2 = [ |/(x)|V(dx)<oo, (A\f 2 }= I Mx)f 2 (x)fi(dx). (1) 

The set X can be taken as the phase space of a particular problem as will 
be the case in this paper. Next we need a finite or infinite orthonormal set 
S = {4> n (x),n = 1, 2, . . . }, selected among the elements of L 2 (X, /i). This set 
spans, by definition, the separable Hilbert subspace Hs and must obey the 
following condition: 

< J\f(x) = \4> n (x)\ 2 < oo almost everywhere . (2) 

n 

Now let us define the family of coherent states {\x)} x( zx in Hs through the 
following linear superposition: 



\x) = -j=—^2<f> n (x)\n) , (3) 

where the states \n) are in one to one correspondence with the functions in 
the set S. This is an injective map X 3 x \— > \x) G Hs (which should be 
continuous with respect to some minimal topology affected to X for which the 
latter is locally compact): These coherent states have two main features: they 
are normalized, (x\x) = 1 and crucially, they resolve the identity in 7ig 

/ \x)(x\Af(x) fi(dx) = I Hs - (4) 
The CS quantization of a classical observable f(x) on X, consists then in 
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associating to f(x) the operator 

/:= / f(x)\x)(x\Af(x)iM(dx). (5) 

This "diagonal" decomposition (in a topological weak sense) may reveal to be 
valid for a wide class of operators. The function f(x) = Af(x) is called upper 
(or contravariant) symbol of the operator / and is non-unique in general. On 
the other hand, the mean value (x\f\x) = Af(x) is called lower (or covariant) 
symbol of /. 



3 Quantization of the particle motion on the circle S* 1 

The motion in the infinite square well potential can be seen as a particular 
case of the motion on the circle S 1 , once we have identified the boundaries of 
the well with each other and imposed Dirichlet conditions on them. Functions 
on this domain will behave as pinched waves on a circle so it is useful to expose 
first the more general case. 

Applying our scheme of quantization we can define the CS's on the circle. The 
measure space X is the cylinder S l x K. = {x = (q,p) | < q < 2tt, p,q e K}, 
i.e. the phase space of a particle moving on the circle, where q and p are 
canonically conjugate variables. We consistently choose the measure on X as 
the usual one, invariant (up to a factor) with respect to canonical transforma- 
tions: fi(dx) = y~ dqdp. The functions (j> n {x) forming the orthonormal system 
needed to construct CS's are suitably weighted Fourier exponentials: 

(t>n{x) = Q V4 e-f ( p -™) 2 e^ n e Z , (6) 

where e > can be arbitrarily small. This parameter includes the Planck con- 
stant together with the physical quantities characterizing the classical motion 
(frequency, mass, etc.). Actually, it represents a regularization. Notice that 
the continuous distribution x i— > |0 n (x)| 2 is the normal law centered at n (for 
the angular momentum variable p) . We establish a one-to-one correspondence 
between the functions <p n and the states \n) which form an orthonormal basis 
of some generic separable Hilbert space 7i that can be viewed or not as a 
subspace of L 2 (X, fj,(dx)). CS's, as vectors in Ti, read then as 
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where the normalization factor 

N{x) = Nip) = \f^T e" e(p - n)2 < oo , (8) 

is a periodic train of normalized Gaussian functions and is proportional to an 
elliptic Theta function. Applying the Poisson summation yields the alternative 
form: 

jV(p) = ]Te 2 ™ p e-^ n2 . (9) 
From this formula it is easy to prove that lim e _ > o Af(p) = 1. 

The CS's dTj) have been previously proposed, however through quite differ- 
ent approaches, by De Bievre-Gonzalez (1992-93) [4|, Kowalski-Rembieliriski- 
Papaloucas (1996) [5], and Gonzalez-Del Olmo (1998) [6]. 



3.1 Quantization of classical observables 



The quantum operator acting on Ti., associated to the classical observable 
f(x), is obtained as in (jHJ). For the most basic one, i.e. the classical observable 
p itself, the procedure yields 

P= N(p)p\p,q)(p,q\iJ l (dx) = ^2n\n)(n\, (10) 



and this is nothing but the angular momentum operator, which reads in an- 
gular position representation (Fourier series): p = — i-^-. 

For an arbitrary function f(q), we have 



/(?)=/ K dx W(p)f(q) \p,q)(p,q\ 



x 

= E e " i(n_n ' )2 Cn-n'(f)\n)(n'\ } (11) 

n,n'& 

where c n (f) is the n-th Fourier coefficient of /. In particular, we have for the 
angular position operator q : 

_±( n _ n ')2 

g = vrl w + i E — (12) 

~r, n — n 



The shift operator is the quantized counterpart of the "Fourier fundamental 
harmonic" : 

= e ~i E \ n + 1 )( n l- ( 13 ) 
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The commutation rule between (|T0|) and (|T3|) gives 



[p,e^] = e l i, (14) 

and is canonical in the sense that it is in exact correspondence with the clas- 
sical Poisson bracket 

{p,e iq } = ie [q . (15) 

Some interesting aspects of other such correspondences are found in [l9j. For 
arbitrary functions of q the commutator 

IP, /5) ] = £(n - n')e-* (n_n ' )2 <*-»'(/) InXn'l, (16) 



can arise interpretational difficulties. In particular, when f(q) = q, i.e. for the 
angle operator 

[p,?]=iE e " i(B_B T \n)(n'\, (17) 

the comparison with the classical bracket {p, q} = 1 is not direct. Actually, 
these difficulties are only apparent if we consider instead the 27r-periodic ex- 
tension to K of f(q). The position observable f(q) = q, originally defined in the 
interval [0, 2tt), acquires then a sawtooth shape and its periodic discontinuities 
are accountable for the discrepancy. In fact the obstacle is circumvented if we 
examine, for instance, the behaviour of the corresponding lower symbols at 
the limit e — > 0. For the angle operator we have 

1 / AT(po-l) ^ e-f^o 

2 V Jv(p ) ; ^ n 



~ 7i + i 



e-*0 n 



where we recognize at the limit the Fourier series of f(q). For the commutator, 
we recover the canonical commutation rule modulo Dirac singularities on the 
lattice 27rZ. 



{po,qo\[p,q]\Po,qo) = l (l + ) 1 ~i + E ^ f n **°* 



£->0 



i + i£%)-27m). (19) 
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4 Quantization of the motion in an infinite well potential 



4-1 The standard quantum context 



Any quantum system trapped inside the infinite square well ^ q ^ L must 
have its wave function equal to zero beyond the boundaries. It is thus natural 
to impose on the wave functions the conditions 

tp(q)=0, q^L and q^O. (20) 

Since the motion takes place only inside the interval [0, L], we may as well 
ignore the rest of the line and replace the constraints (|20|) by the following 
ones: 

^GL 2 ([0,L],dg), #) = ^(i) = 0. (21) 

Moreover, one may consider the periodized well and instead impose the cyclic 
boundary conditions ip(nL) = 0, Vn G Z. 

In either case, stationary states of the trapped particle of mass m are easily 
found from the eigenvalue problem for the Schrodinger operator with Hamil- 
tonian: 

h 2 d 2 

H = H W = -^—. (22) 
2m dar 

This Hamiltonian is self-adjoint [7] on an appropriate dense domain in fl2Tj) . 
Then 

*(q,t) = e-l Ht *(q,0), (23) 
where ty(q,Q) = i/)(q) obeys the eigenvalue equation 

HiP(q) = EiP(q) , (24) 

together with the boundary conditions (I2TI) . Normalized eigenstates and cor- 
responding eigenvalues are then given by 

4>n(q) = y^ sin ( n7r |~) ' ^ q ^ L ' ( 25 ) 
Hij n = E n ij n , n=l,2,..., (26) 



with 



h 2 n 2 2 2 hn 2 2ir 

K= 2^ n S ' UJ= 2^ = Y> (27) 



where T r is the "revival" time to be compared with the purely classical round 
trip time. 
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4-2 The quantum phase space context 



The classical phase space of the motion of the particle is the infinite strip 
X = [0, L] x K = {x = (q,p) | q G [0, L] ,p G M} equipped with the measure: 
p,(dx) = dq dp. A phase trajectory for a given non-zero classical energy = 
\mv 2 is represented in the figure [ZD 

Typically, we have two phases in the periodic particle motion with a given 
energy: one corresponds to positive values of the momentum, p = mv while 
the other one is for negative values, p = —mv. This observation naturally leads 
us to introduce the Hilbert space of two-component complex-valued functions 
(or spinors) square-integrable with respect to p,(dx) : 



Ll 2 {X,»(dx)) ~ C 2 ®L 2 c (X,p(dx)) = U(x) = ( Mx) ), ± e L 2 c {X,n(dx)) 



K(j)_{x)j 



f28) 



We now choose our orthonormal system as formed of the following vector- 



valued functions $ ne (x), K = ±, 



V / \<f> n> -(x)J 
4>n,K{ x ) — v^ce Kp "' ) sin f n7r |^ > k = ± , n = 1, 2, . . . , (29) 



where 



p n = j2mE n = ^-n, (30) 



and the half-width p > is a parameter which has the dimension of a momen- 
tum, say p = find I L with •& > a dimensionless parameter. This parameter 
can be arbitrarily small (like for the classical limit) and, of course, arbitrarily 
large (for a very narrow well, for instance). 

The functions $ njK (x) are continuous, vanish at the boundaries q = and 
q = L of the phase space, and obey the essential finiteness condition ([2]): 



< N{x) = Af(q,p) = AT + (x) + tf_(x) = E E *1,k( x )*»>«( x ) 

ft=± n=l 



-E 

n=l 



OO 

e -^(P-P") 2 +e -^(P+Pn) 2 ,2/ .'/ 



sin 2 ( n7r -jr ) < oo . (31) 
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The expression of jV simplifies to : 

{1 00 1 2 1 

J2 [i-^y^iP-Pn) (32) 
n=— 00 J 

It then becomes apparent that jV and S can be expressed in terms of elliptic 
theta functions. Function S has no physical dimension whereas M has the 
same dimension as c, that is the inverse of an action. 

We are now in measure of defining our vector CS's [8]. We set up a one-to-one 
correspondence between the functions $ n ,«;'s and two-component states 



\n,±) d ^\±)®\n), l+)=( v oJ' l_) = (lJ' (33) 

forming an orthonormal basis of some separable Hilbert space of the form 
/C = C 2 <S> T~L- The latter can be viewed also as the subspace of U^ 2 (X, //(da;)) 
equal to the closure of the linear span of the set of $ nA 's. We choose the 
following set of 2 x 2 diagonal real matrices for our construction of vectorial 
CS's: 

/ 



1\ , , (0 



FJx) 



(34) 



0n,+ (?,P) 

V (j) n ,-(q,p) / 

Note that J\f{x) = E^°=i tr(F n (x) 2 ). Vector CS's, \x,x) <= C 2 <g> H = K, are 
now defined for each 16I and x £ C 2 by the relation 



y.Af(a:) n=l 

In particular, we single out the two orthogonal CS's 

j 00 

|x, k) = ^ F n (x) |n, /c) , « = ± . (36) 

\JjV(x) n=l 

By construction, these states also satisfy the infinite square well boundary 
conditions, namely \x,K) q=Q = \x,n) q= L = 0. Furthermore they fulfill the 
normalizations 

(x, k\x, k) = ^f X ) , J2 (x, k\x, k) = 1 , (37) 



x 



and the resolution of the identity in /C: 



K=± 



I \x)(x\J\f(x)n(dx) = H / / 'F n (q,p)F n/ (q,p)\n,K)(n',K , \dqdp 

JX , , , 1 , , J —00 Jo 



k,k'=± n,n'=l 



\n, k) (n, k\ = a <g) I n = Ifc . (3£ 

k=± n=l 
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where <To denotes the 2x2 identity matrix consistently with the Pauli matrix 
notation to be used in the following. 



4-3 Quantization of classical observables 



The quantization of a generic function f(q,p) on the phase space is given by 
the expression (jSJ), that is for our particular CS choice: 

^ poo pL 

f = J2 / /(<?,p)k,^)(£,«|Af(g,p)dgdp 

~1 J-QO JO 



K=±' 



= \ n )( n '\ 

where 



n,n'=l 



'h 

/-, 



(39) 



f±= d P/ dq(j) n:± (q,p)f(q,p)(p n/t± (q,p) . (40) 
J -oo Jo 

For the particular case in which / is function of p only, f(p), the operator is 



given by 



^ poo rL, 

f = / / f(j>)\x,K)(x,K\Af(q,p)dqdp 

K=± J-oo JO 

( 



■y 00 

— ?= \ n ){ n \ ® 



/+ 



(41) 



with 



/± = J dpf(p) exp ( - — (p =f p„) 2 ) . (42) 



p- 

Note that this operator is diagonal on the |n, «) basis. 



^.5.1 Momentum and Energy 

In particular, using /(p) = p, one gets the operator 

p = ® |n)(n| , (43) 

n=l 



where (73 = ( -1 ) * s a Pauli matrix. 



For /(p) = p 2 , which is proportional to the Hamiltonian, the quantum coun- 
terpart reads as 

— o 2 00 o 2 

= + E^o ® |n)(n| = ^-I K + (p) 2 . (44) 

Z n=l A 
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Note that this implies that the operator for the square of momentum does not 
coincide with the square of the momentum operator. Actually they coincide 
up to 0(h 2 ). 



4.3.2 Position 

For a general function of the position f(q) our quantization procedure yields 
the following operator: 



1 2 

/= £ e"^-^ [<*„-«'(/) ®\n)(n'\ 



n,n'=l 



where 



In particular, for /(g) = q we get the "position" operator 



(45) 



(46) 



_ L 2L ~ _ 1 (Pn _p ,)2 
q = -Ik 9" V e 

2 n,^l, 

n + n'=2fc + l 



(n — n') 2 (n + n') 2 



(47) 

with fcGN. Note the appearance of the classical mean value for the position 
on the diagonal. 



4-3.3 Commutation rules 

Now, in order to see to what extent these momentum and position opera- 
tors differ from their classical (canonical) counterparts, let us consider their 
commutator: 



2h 00 

[q,P} = — Y] Cn,n'0- 3 ®\n){n'\ 

K -A I 

n+n' = 2fc + l 



e 4 p- 



p 2 n (n — n ) 



(n — n') 2 (n + n') 2 



(48) 
(49) 



This is an infinite antisymmetric real matrix. The respective spectra of finite 
matrix approximations of this operator and of position and momentum op- 
erators are compared in figures [2] and [3] for various values of the regulator 
p = hnfi/L — 1? in units h = 1, L — 7r. When p takes large values, one can 
see that the eigenvalues of [q,p] accumulate around ±i, i.e. they become al- 
most canonical. Conversely, when p — > all eigenvalues become null, which 
corresponds to the classical limit. 
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4-3.4 Evolution operator 



The Hamiltonian of a spinless particle trapped inside the well is simply H = 
p 2 /2m. Its quantum counterpart therefore is H = p 2 /2m. The unitary evolu- 
tion operator, as usual, is given by 



U(t) = e~* m = e'™** e'^o-Q <g> \n) (n\ . (50) 

71=1 



Note the appearance of the global time-dependent phase factor with frequency 
ud which can be compared with the revival frequency 

UJ » = ^L 2= —- (51) 



5 Quantum behaviour through lower symbols 

Lower symbols are computed with normalized CS's. The latter are denoted as 
follows 

\x) = \x, +) + \x, -) . (52) 
Hence, the lower symbol of a quantum observable A should be computed as 

A(x) = (x\A\x) = A ++ (x) + A + _(x) + A^ + (x) + A__(x) . 

This gives the following results for the observables previously considered : 

5.0.5 Position 

In the same way, the mean value of the position operator in a vector CS \x) 
is given by: 

( x \q\ x ) = \ ~ Q(Q,P) » ( 53 ) 

where we can distinguish the classical mean value for the position corrected 
by the function 



n,n =l,n?=n 
n+n'=2k+l 



1 



(n — n') 2 (n + n'y 



x 



x 



e -^[(P-Pn) 2 +(P-JV) 2 ] + 



_|_ e -^[(P+P«) 2 + (P+Pn') 2 ] 



sin ( nn j^j sm ( n ' 7T J^) ' 
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This function depends on the parameter d as we show in figure O with a 
numerical approximation using finite matrices. As for p, we calculate the dis- 
persion defined as 

AQ = ^q 2 - q 2 . (55) 
Its behaviour for different values of d is shown in figure Q3 



5.0.6 Time evolution of position 

The change through time of the position operator is given by the transforma- 
tion q(t) := U^{t) qU(t), and differs from q by the insertion of an oscillating 
term in the series. Its lower symbol is given by 



(x\q(t)\x) = - - Q(q,p,t) , 



(56) 



where this time the series have the form 



Q{q,p,t) 



2L 1 



E 



e -2^(p'-^)* e "3^ (Pn " p "' )2 x 



n,n'=l,n^n' 
n+n'=2k+l 



X 



sin nir 



sin n ix— 



(n — n') 2 (n + n') 2 



X 



(57) 



Note that the time dependence manifests itself in the form of a Fourier series 
of with frequencies (n 2 — n' 2 ) fm 2 /2mL 2 . This corresponds to the circulation 
of the wave packet inside the well. 



5. 0. 7 Momentum 

The mean value of the momentum operator in a vector CS \x) is given by the 
afline combination: 



(x\p\x) 



M(x) 



X 



M{X) =cY,Pn 
n=l 



2 / 9 

sin j rnr — 



(5f 



This function reproduce the profile of the function p, as can be seen in the 
figure [5j We calculate then the dispersion AP, defined as 



AP 



p2 _ p2^ 



(59) 



using the mean values in a CS \x). Its behaviour as a function of x is shown 
in figure [71 
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5.0.8 Position-momentum commutator 



The mean value of the commutator in a normalized state ^ = y ^_ J is the 
pure imaginary expression: 



2ih 



oo 




TV 



n — n 



x 



n + n> ' = 2fc + l 



1 



1 



S((0 + |n)(n> + )-(0_|n)(n'|0-))- (60) 



x 




(n + n') 2 



Given the symmetry and the real-valuedness of states (1361) . the mean value of 
the commutator when \l/ is one of our CS's vanish, even if the operator does 
not. This result is due to the symmetric spectrum of the commutator around 
0. As is shown in Part c) of figures [21 the eigenvalues of the commutator 
tend to ±ih as p, i.e. increases. Still, there are some points with modulus 
less than h. This leads to dispersions AQAP in CS's \x) that are no longer 
bounded from below by h/2. Actually, the lower bound of this product, for 
a region in the phase space as large as we wish, decreases as $ diminish. A 
numerical approximation is shown in figure El 



6 Discussion 

From the mean values of the operators obtained here, we verify that our CS 
quantization gives well-behaved momentum and position operators. The clas- 
sical limit is reached once the appropriate limit for the parameter •& is found. If 
we consider the behaviour of the observables as a function of the dimensionless 
quantity $ = pL/hn, at the limit i? — > and when the Gaussian functions for 
the momentum become very narrow, the lower symbol of the position opera- 
tor is q ~ L/2. This corresponds to the classical average value position in the 
well. On the other hand, at the limit d — > oo, for which the involved Gaussians 
spread to constant functions, the mean value {x\q\x) converges numerically to 
the function q. In other words, our position operator yields a fair quantitative 
description for the quantum localization within the well. The lower symbol 
(x\p\x) behaves as a stair-step function for p close to and progressively fits 
the function p when p increases. These behaviours are well illustrated in the 
figures [H and [51 The effect of the parameter $ is also noticeable in the dis- 
persions of q and p. Here, the variations of the full width at half maximum of 
the Gaussian function reveal different dispersions for the operators. Clearly, if 
a classical behaviour is sought, the values of $ have to be chosen near 0. This 
gives localized values for the observables. The numerical explorations shown 
in figures [6] and [7] give a good account of this modulation. Consistently with 
the previous results, the behaviour of the product AQAP at low values of $ 
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shows uncorrelated observables at any point in the phase space, whereas at 
large values of this parameter the product is constant and almost equal to the 
canonical quantum lowest limit h/2. This is shown in figure [751 

It is interesting to note that if we replace the Gaussian distribution, used here 
for the p variable in the construction of the CS's, by any positive even proba- 
bility distribution K. e p \— > w(p) such that J2 n w i.P ~~ n ) < 00 the results are 
not so different! The momentum spectrum is still Z and the energy spectrum 
has the form {n 2 + constant}. In this regard, an interesting approach com- 
bining mathematical statistics concepts and group theoretical constructions 
of CS's has been recently developed by Heller and Wang pZOfiZTj . 

The work presented here has possible applications to those particular physical 
problems where the square well is used as a model for impenetrable barriers 
|22j . in the spirit of what has been done in [23] . 

The generalization to higher- dimensional infinite potential wells is more or less 
tractable, depending on the geometry of the barriers. This includes quantum 
dots and other quantum traps. Nevertheless, we believe that the simplicity 
and the universality of the method proposed in the present work should reveal 
itself useful for this purpose. 

Author Garcia de Leon wishes to acknowledge the Consejo Nacional de Ciencia 
y Tecnologia (CONACyT) for its support. 
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Fig. .1. Phase trajectory of the particle in the infinite square- well. 
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Fig. .2. Eigenvalues of q, p and [q, p] for increasing values of the characteristic 
momentum p = hirfi/L of the system, and computed for N x N approximation 
matrices. Units have been chosen such that H = 1, L = it so that p = ~Q and p n = n. 
Note that for q with p small, the eigenvalues adjust to the classical mean value L/2. 
The spectrum of p is independent of p as is shown in (03]). For the commutator, the 
values are purely imaginary. 
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Fig. .3. Continued from figure [?2l N x N approximation matrices eigenvalues of q, 
p and [q,p] for increasingly larger values of p = h-n-d/L = $ in units h = 1, L = tt. 
The spectrum of p is independent of p as is shown in (03]). For the commutator, the 
eigenvalues are purely imaginary and tend to accumulate around ih and —ih as p 
increases. 
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Fig. .4. The lower symbol q depicted for various values of p = hwd/L = i? in units 
H = 1, L = it. Note the way the mean value fits the function q when p is large, and 
approaches the classical average in the well for low values of the parameter. 
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Fig. .5. The lower symbol p depicted for various values of p = hir&/L = -d in units 
H = 1, L = 7T. The function becomes smoother when p is large. 
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Fig. .6. Variance of q depicted for various values of p = hnd/L = d in units h = 1, 
L = 7r. Note how different dispersions are revealed just by changing the width of the 
Gaussian function of the p variable. Low dispersion, close to classical, is found for 
$ near and the quantum behaviour is recovered at large values of the parameter. 
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Fig. .7. Variance of p depicted for various values of p = hn-d/L = i? in units H = 1, 
L = 7r. Consistently with g, a well localized momentum is found for low values of 
the parameter. This is actually expected since the Gaussian becomes very narrow. 
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Fig. .8. Product AQAP for various values of p = fm-djL = i9 in units h = 1, L = ir. 
Note the modification of the vertical scale from one picture to another. Again, the 
pair position-momentum tends to decorrelate at low values of the parameter, like 
they should do in the classical limit. On the other hand it approaches the usual 
quantum-conjugate pair at high values of p. 
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